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ABSTRACT. The Friedman universe is re-examined in a context that is non-standard only in 
that the properties of matter are postulated in the form of an action principle. Applications 
to equilibrium configurations of ideal stars have already been reported. In this paper we apply 
the same theory to a fresh examination of the Friedman universe. The results agree with 
standard theory in the case of low densities. A suggestion is made to replace "vacuum energy" 
by "external force" . 

1. Introduction 

To make this paper reasonably self contained, we present a brief summary of the main 
idea and its merits, though it was done to the best of our ability in [F] . 

Very briefly, this work is driven by a concern about the standard approach to rela- 
tivistic dynamics that is presented in the textbooks, as a fusion of General Relativity and 
Nonrelativistic Thermodynamics. The basis for it all is Tolman's formula for the energy 
momentum tensor. 



used as the source of Einstein's field equations. It involves 2 scalar fields, and a vector 
field U that is dismissed by an argument that involves "co-moving coordinates", and a 
normalization condition that we believe had better been replaced by an equation of motion. 
There is also a baryon density that seems curiously unrelated to the rest. 

Our approach takes General Relativity with Einsteins equations for granted, and tries 
to find model sources that have reasonable properties. Classical thermodynamics is not 
posited as input, instead the interpretation is to be extracted from the model itself. The 
fact that this interpretation can be understood in thermodynamcal terms results from an 
examination of the structure of the model, it is not postulated independently. 
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The four degrees of freedom that are represented by the velocity field in the standard 
treatment are not fully exploited in applications to stellar dynamics or cosmology. The 
space components are removed by specialization to static phenomena, characterized by the 
absence of flow. An alternative approach avoids introducing an independent velocity held 
by focusing the attention on irrotational flows, as flrst proposed by Jeans in 1902. The 
possibility of describing turbulence is thus discarded from the beginning; this is a strategy 
that is standard in non-relativistic hydrodynamics whenever laminar flow is expected to 
dominate, as in the theory of propagation of sound or heat. Thus in nonrelativistic hydro- 
dynamics matter is described in terms of a density and a velocity potential. The equations 
supposedly satisfled by these flelds are Bernoulli's equation (derived from Newton's equa- 
tion) and the continuity equation. The theory admits a formulation in terms of an action 
principle. The proposed relativistic version of this theory involves a scalar held p and a 
velocity potential '0, 

The nonrelativistic theory is obtained by setting 

and identifying $ with the ordinary velocity potential. The equations of motion are 

^(^'^^V',mV',.-c^) = ^ (1-2) 

and the equation of continuity, 

d^iV^pg^^''iP,.) = 0. (1.3) 

It is natural to interpret V as an internal energy (including the rest mass), this leads to 
an identification of the lagrangian density with pressure 

dp 

Back to the standard theory. After disposing of the space components of the velocity 
field it is the turn of the time component. This embarrassing quantity is always eliminated, 
and always without comment, by fiat: gi_ii,U^U^ = 1! In the absence of an action principle 
this does not lead to internal contradictions; our discomfort arises from the fact that there 
is no justification and every reason to be worried. The original context of this condition is 
the motion of a particle along a geodesic. Here — dx^ /ds and g^vU^U^ is a constant 
of the motion. Assigning the value unity is the relativistic version of the identification 
E = H in classical mechanics. But the condition is not legitimate when more than one 
particle is involved, unless each particle moves independently of the others; that is, unless 
there is no interaction other than gravity. If other forces are present it is absurd to impose 
this condition on each particle separately. In fact, in relativistic particle theory one does 
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not know what to do. The Bethe-Salpeter equation is an unsatisfactory way to deal with 

the problem, for strictly it makes sense only in the static approximation. Rclativistic 
mechanics of interacting particles does not exist; the answer is field theory. Perhaps the 
best that can be said for the condition g,j,vU^U^ = 1 is that, in the static case, when it 
reduces to goQ{U^Y = 1, it is really irrelevant because this factor can be absorbed into 
the density. True, but this merely transfers the problem to the density; which density is 
expected to satisfy an equation of state? The alternative approach replaces by g'^'^t/j^i, 
and Eq.(1.2) constrains the invariant. In the static case the derivative of this equation 
with respect to the radial coordinate is exactly the same as the hydrostatic condition of 
the standard approach. The only new element here is a new boundary condition that turns 
out to be very beneficial in stellar dynamics. [F] 

The other equation, Eq.(1.3), is a conservation law. In the standard theory there is 
no conserved quantity. The density in Eq.(l.l) is not conserved and its interpretation is 
controversial. For example, an equation of state relating p to p, proposed by Chandrasekhar 
[Cl], was taken over by other authors [OV] who replaced p by the energy density p — p. 
The framework is sometimes expanded with the addition of a conserved baryon density; it 
does not have natural place in the theory, but its mere existence was used to good effect 
to justify the introduction of an equation of state [C2]. 

Nevertheless, the conservation law (1.3) does not impose additional constraints on the 
theory, for when Eq.(1.2) holds it is equivalent to the contracted Bianchi identity. 

From now on, in this paper we are using units such that c = 1. 

Summary and conclusions 

In Section 2 we derive the equations that govern the metric and matter fields in 
the Friedman universe. The metric is assumed to be such that, in a suitable coordinate 
system, it takes the form of the Robertson- Walker metric. As in the standard theory we 
seek solutions without fiow; that is, if) depends only on the (Robertson- Walker) time. The 
equations that result from this, for the metric and for the fields p,p and ip, would agree 
with the corresponding equations of the standard theory if one would replace the equation 

hy ip = 1 (and replace p hy p + p). Of course, the model does not allow that, but it is a 
good approximation for very low pressure and weak gravitational fields. 

The time derivative of this equation (1.4) is a consequence of the contracted Bianchi 
idcnties and appears in the traditional approach as well. Perhaps the most important 
property of the model is that Eq.(1.4) is stronger than its time derivative; it fixes an 
additional constant. It can be used to determine the value of ip at the present time. This 
will be done in Section 3 after selecting an equation of state. 

In Section 3 we also present results of numerical integration of the equations, forwards 
and backwards in time from present conditions, using a polytropic equation of state. 

The conclusion is that there is very little difference between the predictions of the 
respective theories. The merit of our proposal, seen in the present limited context, is, in 
our opinion, aesthetic and improved logical coherence. Most remarkable is the insensitivity 
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of the results with respect to wide variation of the polytropic index. This circumstance 
gives some confidence that the results are not overly sensitive to the choice of the equation 
of state. 

In Section 4 we present a suggestion for overcoming the apparent discrepancy between 
the observed density of matter of the universe and the critical density favored by the 
measured value of the Hubble constant and other recent data. 

More 

To include non-laminar flows, a generalization is needed, leading to the diflicult prob- 
lems including turbulence. Accommodating rotating stars and the Kerr metric is an easier 
task and it may serve as a first step towards a general theory. Quantization of the model, if 
it proves possible, may throw some light on fluctuations. And thinking about this problem 
might possibly lead to some inspiration about quantum gravity. 



2. The equations 

Einstein's equations, 

GfjLv + Afif^j^ = SttGTh^, 

will be applied to the following situation. 

It is assumed that there are coordinates t, r, 6, (p in which the line element takes the 

form 

ds^ = dt" - R\t){^^ + r'^dn'). 

The components of the Einstein tensor depend only on so the same must be true of the 
energy momentum tensor. The matter action is 

^matter = ^ (^'^^.^V,. - 1) - ^[p]) =: / d^x^ C. 

and the associated energy momentum tensor 

T^,u = Pi^,^ti^,jy - g^iu^- (2.1) 
The flelds p, ij) are subject to the Euler-Lagrange equations 

d^{^pg^^^,.) = 0, \{g^''^,^i^,. - 1) = ^ • 

We shall assume that both flelds depend on t only, then these equations reduce to 

|(«3^)=0, l#-l) = f. (2.2) 

The components of the Einstein tensors may be found in the textboks [W]; the nonzero 
components are 

Gtt = {3/R^){R^ + k), Gij = -R-\2RR + R^ + k)gij, 
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where R^g is the spatial part of the metric. Einstein's equations boil down to 

(3/i?2)(i?2 + A = 877G{piP^ -p), 

and one of the Bianchi identities (equation of 'energy conservation'). But the latter is 
implied by the wave equations and can be ignored. 

The energy momentum density T^,^ = y/^T/j,,^ obtained from Eq.(2.1)(2.1) satisfies 

(t;),, + V^,^c + l^/^pg'^fM,f3 = 0, fi = t,r, e, 4>. 



2 

In the present circumstances this reduces to one equation, 

R^p=^^{R'pi^% (2.4) 

If we replace t/^ by 1 in Eq.s (2.3) and (2.4), and change the notation by translating 
p ^ p, ip'^p p + p we obtain the equations that define the Friedman universe according 
to the standard theory. 

One can also derive the last result directly from the wave equations (2.2). The first 
allows to define a constant k, 

R^ptp = K, 

and the second gives upon derivation with respect to the time the result 

p/p = (2.5) 

hence R^p = R^pip^/^ ^ Kti: ^ {d/dt){R^pip'^). 
To summarize, the full set of equations is 

R-\R^ + k) = ^{p^'-p) + A/3, R'p^ = K, ^(^^-1) = ^ (2.5-7) 

to be completed by a choice of the functional V[p]; that is, an equation of state. Note that 
(2.7) is stronger than (2.5) 

In the traditional approach the equations are instead 

R-' {R' + k) = ^p + A/3, R'p =l{R\p + p))=0. (2.8-9) 

Our result differs from the standard theory in several minor respects. To begin with, 
Eq.(2.7) fixes ip; that here plays the role of U^. It is not unity, but it does not differ 
significantly from unity if the pressure is negligible. If we replace this equation by its time 
derivative, then as was seen we end up with the standard set of equations by regarding 
the energy density pip'^ — p as the fundamental variable and identifying it with the field p 
of the standard theory. The density used in the equations of state are not the same, but 
that is unimportant if the pressure is negligible. 
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We pointed out that the Bianchi identity (2.4) can be derived from the Euler-Lagrange 

equations (2.2) by differentiating the second one with respect to to the time. The more 
important difference between our model and the standard approach is that the equation 

- 1) = f (2.10) 

contains information that is lost upon differentiation. To understand what is involved here 
let us consider the case of of a polytropic equation of state, V = ap'^, with 7 = 1 + 1/n. 
Then this equation reduces to 

^ = (2.11) 
p 2n + 2 ^ ^ 

Our model predicts the present value of ip in terms of the density and pressure. Unfor- 
tunately, the value exceeds unity by an amount that is far too small to bridge the gap 
between observed density and critical density. Some more discussion of this point may be 
found in Section 4. 

Finally, the appearance of dV/dp in the equations has the effect that, though we can 
eliminate the field ip, it is not enough to make a guess about the probable dependence of 
pressure and density on R{t). So we are in a somewhat more constrained situation in that 
we shall need an equation of state. 

We shall explore the possibility that the universe is a polytrope; that is, there is a 
constant 

7 = 1 + - 

n 

such that 

V[p] = ap'^ , p = {a/n)p'^ . 
Then there is a function / (the Emden function) such that 

V = ar+\ p = r, p={a/n)fp. 

The value of / is closely related to the temperature. Eq.s (2.8-9) become 

R-\ie + k)=^-^p{^^-^)) + K/2,, R'p^^K, V''-l = 27a/. 
6 n 

Let po,Ro,... denote the values of the variables at the present time; we shall take 
i?o = 1 and also ^0 = I5 so that the unit of time is the Hubble time Rq/Rq. Then 

1 = -i^PoiA ) + A - /c, 

where A is the reduced cosmo logical constant. Since the pressure pq is neglible we have to 
a good approximation ipQ — afo/n — 1. Set p — ppcr-, then 

k = poiipo - fo/n) + X-1 
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and our equations take the form that will be used for computations, 



l^iip^ - af/n) + A -k, R^fxip = //o^o, - 1 = 2^af. (2.12 - 14) 



Since it is difficult to estimate the value of the parameter a we shall set 

afit) = bg{t), go = 1, ii/iio = 

and consider the constant 6 as a free parameter. 

We shall compare the results with the standard model, where the equations are 



i?2 = + - k, -R^iiaf = ^ {R^fi{l + af/n)) 



d 



n ' " dt 

Neglecting the effect of pressure gives complete agreement. 



3. Numerical integration 

We have already set c = 1. As unit of time we have taken Hubble time, Rq/Rq. which 
is of the order of the current estimate of the age of the universe, 18 billion years. The 
parameters are n, |Uo, A, h. At time t = the pressure is relatively unimportant, and the 
function ip is close to unity. We begin our investigation with the case //o = l,A = 0, h = 
10-^, initial conditions a/(0) = 6, R{Q) = 1. 

For small values of a, when the pressure is negligible at most times, there is no 
difference between the standard theory and the model. Since the pressure is negligible at 
present, and diminishes in the future, there is no significant difference between the two 
theories in the future. The function R{t), renormalized is 1 at t = (present time) and 
vanishes at time —Iq (big bang). 

The most interesting value of the polytropic index is n = 3, which is appropriate for a 
matter dominated universe and hence for most of the past and all of the future. We start 
from the reference point 

(n,//o,A,6)-(3,l,0,10-«) 
and vary one parameter at a time. At the refence point we have 

to = -.666, qo = -R{t) = .5. 

We shall assume that experiment favors a greater age and a smaller value of go- 

Increasing n has negligible effect, but smaller values of n makes the universe younger. 
Thus for n = .2, go remains at .5 and 

to = —.365, (standard theory), to = —.596 (dynamical model). 

Increasing the pressure has little effect unless it is very high. Thus for b — .1, 

to = —.611, go = .56 (st. th.), to = —.575, qo = .67 (dyn.mod) 
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This is going in the wrong direction, but probably a reasonable mean value for the param- 
eter b is much less than .1. 

By current experimental evidence the value of no is much less than unity. We consider 
the alternative reference point 

(n,//o,A,6) = (3,.3,0,10-^). 

Then 

to = —.809, qo — .15 (either theory). 
After increasing the pressure as before, with 6 = .1, 

to = — -741, = -2 (either theory). 

From the experimental point of view there is very little to distinguish between the two 
theories. 

4. An alternative to vacuum energy 

Recent observations point to a spatially flat universe (/c = 0), while direct observation 
of the matter density gives only 30 percent of critical. The only popular way to resolve the 
discrepancy is to attribute the missing 70% to "vacuum energy" represented by a positive 
cosmological constant. A theoretical prejudice against a positive cosmological constant 
(no invariant vacuum and consequently no supersymmetry) is a strong incentive to look 
for alternatives. 

The factor ip'^ that accompanies the density p in Eq.(2.5) seems to offer an opportunity, 
but it is far too close to unity to be helpful. A generalization of the model may be the 
answer. 

Scale transformations in the traditional approach 

We have alluded to the fact that fixing the normalization of the four-velocity field U 
amounts to fixing the interpretation of the density. Here we shall be more precise. 
Let qhe a, real, positive constant and define U,p,p by 

U = qU, p = q~^p, p = p. 

Tolmans formula goes over to T^jy = pUnUu — pQnu', it is scale invariant. The scale is, 
nevertheless, meaningful. This for various obvious physical reasons and also since the 
transformation affects the equation of state. For simplicity, suppose it has the form p = 
ap"'] then p = ap'^ , with d = q~'^'^a. If the value of the constant a is known, theoretically 
or experimentally, then it implies a definition of density that fixes the scale. If we had no 
such knowledge then the density would not be defined and setting either g^'^U^Ui, = 1 or 
g^'^U^Uv = 1 would be equally justified. 
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A cosmic force 

What would be the effect of adding a term proportional to the density to our La- 
grangian? We take 

>C = I (^'^^VmV'. -l)-V-Fp. 

In the nonrclativistic theory the addition of the term proportional to the density is not 
observable, for it corresponds to a sort of gauge transformation ([FW], page 301). If 
•0 = c^t + $, then g^^ip"^ — 2c^$. Here $ is the non-relativistic velocity potential; 
adding a constant to it does not affect the velocity y^- 

The extra term can be eliminated by a scale transformation 

Our Lagrangian is not invariant; it takes the form 

>C = f {g^'-i^f.i^. - q-^) -V- {F/q^)p. 
If we choose q such that F = (g^ — l)/2 then it becomes 

>C = I [g^'-i^^.i^. -l)-V. 

In the cosmological context the pressure is ill defined and usually neglected in the case of 
a matter dominated universe. Consequently, there is no independent determination of the 
external force or, what is the same, the scale of the density. The wave equations now give, 

in the case of a neglible pressure, •0 pa 1. The effective source for the gravitational field is 
p = p/q^. If p is identified with the mass density as measured locally, then a value = 3, 
or F = 1, will appear to solve the problem of the missing mass/energy and there will be 
no need to invoke a non-vanishing cosmological constant. 

Calling in a cosmic force is not a complete solution to the problem, since one would 
need to explain why it does not play a role in the local context. The inverse is perhaps 
more plausible, a term proportional to the density may manifest itself locally, leading to a 
low estimate of the mass density. The factor q^ would then have the meaning of a "biasing 
parameter" [SW]. 
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